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ABSTRACT 



Context. The Milky Way bulge is the nearest galactic bulge and the best laboratory for studies of stellar populations in spheroids based 
on individual stellar abundances and kinematics. These studies are challenged by the strongly variable and often large extinction on a 
small spatial scale. 

Aims. We use the Vista Variables in the Via Lactea (VVV) ESO public survey data to measure extinction values in the complete area 
of the Galactic bulge covered by the survey at high resolution. 

Methods. We derive reddening values using the method described in Paper I. This is based on measuring the mean (J - K s ) color of 
red clump giants in small subfields of 2' x 2' to 6' x 6' in the following bulge area: -10.3° < b < +5.1° and -10.0° < I < +10.4°. To 
determine the reddening values E(J - K s ) for each region we measure the RC color and compare it to the (/ - K s ) color of RC stars 
measured in Baade's Window, for which we adopt E(B - V) = 0.55. This allows us to construct a reddening map sensitive to small 
scale variations minimizing the problems arising from differential extinction. 

Results. The significant reddening variations are clearly observed on spatial scales as small as 2'. We find a good agreement between 
our extinction measurements and Schlegel maps in the outer bulge, but, as already stated in the literature the Schlegel maps are not 
reliable for regions within \b\ Sj 6°. In the inner regions we compare our results with maps derived from DENIS and Spitzer surveys. 
While we find good agreement with other studies in the corresponding overlapping regions, our extinction map has better quality due 
to both higher resolution and a more complete spatial coverage in the Bulge. We investigate the importance of differential reddening 
and demonstrate the need for high resolution extinction maps for detailed studies of Bulge stellar populations and structure. 
Conclusions. We present the first extinction map covering uniformly ~ 315 sq. deg. of the Milky Way bulge at high spatial resolution. 
We show that, in a 30 arcmin window, at a latitude of b = -4°, which corresponds to often studied low extinction window, the so 
called Baade's Window, A Ks values can vary by up to 0.1 mag. Larger extinction variations are observed at lower Galactic latitudes. 
The extinction variations on scales of up to 2'-6' must be taken into account when analysing the stellar populations of the Bulge. 

Key words. Galaxy: bulge - ISM: dust, extinction - stars: abundances 
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ctf 1. Introduction 

The Bulge is one of the major components of the Galaxy and 
it can be studied at a unique level of detail, in comparison to 
those of external galaxies, thanks to its proximity which en- 
ables properties of individual stars to be measured. Bulge stars 
hold the imprints of how our Galaxy formed and evolved. The 
Bulge star formation and chemical enrichment history can be 
investigated in d etail through analysis of stellar ages and chem 
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Bensb v et al 
l2011bl) . 

Our current understanding of the Milky Way Bulge struc- 
ture and stellar content can be summarized as follows: i) 
the B ulge structure is dominated by the Bar (IStanek et al.l 
1994) which appears to have a peanut or X-shape in the 
outer regions dMc William & Zoccalil l2010t iNataf etal.1 l2010t 
ISaito et aHl201 lh: ii) there is a radial stellar metallicity gradi- 
ent (IZoccali et al.l 120081: iJohnson et ail 1201 ll) th at most likely 
flattens in the inner regions (IRich et al.l 120071) : iii) in addi- 
tion to the population associated with the bar, there are hints 
for a second chemically and kinematically distinct component 
(Babusiau xltafll201rt iHill et alJl201 lb IGonzalez etafll2011bt 
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jens bv et al. 1201 lh: and iv) the bulge is dominated by old star s 
( ZoccalietalJ 120031 iBrown et al l lioidt Iciarkson et al.1 1201 lh . 
although there are indications for a possible younger popula- 
tion of metal-rich st ars as traced by Bulge microlensed dwarfs 
dBensbv et alj|201 lh . However, the general view of these prop- 
erties, based on a more complete coverage, is expected to pro- 
vide the final picture of the Milky Way bulge formation his- 
tory. Several projects aiming to obtain this general view of 
the Bulge are either ongoing or planned for the near future 
(e.g.lMinniti et alJl20ld lFeltzind[20TTt iNess & Freemanll2012t 
Elia Garc ia Perez et al.ll2012l) . The interpretation of these obser- 
vations is challenged by the variations in the dust extinction 
properties on small scales towards the different Bulge regions. 

The currently available extinction maps of the Bulge regions 
have different resolutions and coverag es. The most commonl y 
used extinction maps were those of Burstein & Hei lesl (Il982l) . 
later superseded by Schlegel maps (ISchlegel et al. 1998). The 
Schlegel extinction maps are full-sky maps of the dust color 
temperature based on IRAS and DIRBE experiments which are 
then normalized to E(B-V) values using a calibration of colors 
of background galaxies. Unfortunately, these maps suffer from 
large uncertainties in regions towards the Bulge. Specifically, as 
stated in the appendix C of Schle gel et al.l d 1998b . the tempera- 
ture structure of the Galaxy towards latitudes \b\ < 5° is not well 
defined and contaminating sources have not been completely re- 
moved. Furthermore, as we show below, the resolution of these 
maps is too coarse. As a result, the Schlegel reddening values 
towards the Bulge are unreliable. 

Clearly, an extinction map covering \b\ < 5°, with sufficient 
resolution to resolve the small spatial scale variations of extinc- 
tion, is important for an understanding of the general properties 
of the inner Galactic regions. Using red giant bra nch (RGB) stars 
from t he near-IR photometric survey DENIS, ISchultheis et al.l 
( 1999) provided a high resolution (2') map in the inner regions 
of the Bulge (\b\ < 2°). This map provides the sufficient resolu- 
tion to analyse such inner Galactic regions, but its coverage is 
too limited to permit the study of the global Bulge stellar popu- 
lations and structural prop erties. 

More recently, Marshal l et al.l (120061) provided a full 3D ex- 
tinction map of the Bulge. 3D maps are indeed of high impor- 
tance for Galactic studies towards the inner Galaxy as they hold 
the additional distance information. However, in order to build 
3D extinction maps, a sufficient number of stars is required in 
each resolution element, which implies that only a modest re so- 
lution can be achieved. The method in Marshal l et al.l (120061) is 
based on a comparison of 2MASS photometry to the Besancon 
model. At low Galactic latitudes (\b\ < 3°) 2MASS suffers 
strongly from incompleteness and blending due to too low res- 
olution and high stellar density, limiting the ability to derive 
correct star counts and extinction. The resolution of Marshall's 
3D extinction map is 15', which is again too coarse to fully 
overcome the differential reddening problems at lower latitudes 
\b\ < 4°. 

Other important extinction maps are available with a par- 
tial c overage of t he Bulge such as those of iKunder et al.l (2008) 
and ISumil (12004 . However, for detailed studies of the Bulge, 
a high resolution and homogeneous extinction map covering 
the complete area of interest is missing. The deep and near- 
IR photometry of t he VVV ESO public survey of the Bulge 



IK pnotometry or tne VVV tsu public survey or tne tsulge 
dMinniti et al. 20101) provides t he ideal dataset to create such a 
map. In|Gonzalez et al. (201 lc, Paper I) we showed an effective 
technique based on the color of red clump (RC) stars to trace 
the effects of extinction. In this article, we extend our analysis 
to the complete Bulge region covered by the survey and derive 



the first complete extinction map of the Bulge, using a homoge- 
neous population, namely the red clump, at a resolution of 2'-6' 
(Sect.0. We compare our map to extinction maps in the litera- 
ture in the corresponding overlap regions (Sect. 2), and we dis- 
cuss the importance of high resolution extinction map in future 
Bulge studies (Sect. [5]). 



2. Determination of A Ks values 

2.1. VVV survey data 

In this work we use the ~ 315 sq. deg. coverage of the Galactic 
bulge from the VVV survey. This correspon ds to the first data re- 
lease (DR1) of the survey, as described in lSaito et al.l (1201 2b . In 
Paper I, we presented a detailed description of a method used to 
produce the multiband (J, H, and K s ) catalogs required for our 
global study of the Bulge stellar populations. Thus, we now pro- 
vide only a brief summary of the adopted procedure. Single band 
catalogs for each VVV image, so-called tiles, are produced at the 
Cambridge Astronomical Survey Unit (CASU). These catalogs 
are then individually calibrated to the 2MASS photometric sys- 
tem by comparison of high quality stellar detections from both 
catalogs, in the magnitude range between 12 < K„ < 13 mag 
where both catalogues overlap. Sources brighter than K s = 12 
mag in each VVV tile are removed and are then replaced with 
2MASS stellar detections, that have been flagged with high qual- 
ity photometry. This ensures the required photometric agreement 
between the 2MASS and VVV surveys and the correction of sat- 
uration at the bright end (K s < 12) of the VVV catalogs. These 
single band catalogs are then matched using STILTS (Tavlor 
2006) to produce our final multiband catalogs. 

2.2. Mean color of the Bulge red clump 

Studies of the chemical abundances and age measurements of the 
Bulge towards Baade's Wind ow, have shown tha t the Bulge is 
predominantly old (~10 Gvr) dOrtolani et al.l 1 9951: IZoccali et all 
120031: IBrown et al1l2010t Iciarkson et al.ll201 ll) and that it shows 
a broad metall icity distribution with a peak near Solar metallic- 
ity dMc William & Richlll994tlZoccali et alJl2008l:lJohnson et ail 
201 1). Stellar evolutionary models can thus provide specific con- 
straints on the expected intrinsic color of stars along the different 
evolutionary stages. The observed color can then be compared to 
that expected from models and its difference can be linked to an 
external factor, such as extinction. The metallicity of the Bulge 
is not spa tially uniform, but it shows a gradient along its mi- 
nor axis dZoccali et al. 2008). Furthermore, the extent of these 
gradients towards other regions is not yet clear. For this reason, 
uncertainties may arise from direct comparison of the observed 
red giant branch (RGB) stars' colors with respect to those ex- 
pected from models. Fortunately, this effect can be minimized 
by the use of RC giants, as the mean color of these stars has a 
small dependance on these parameters. 

In Paper 1 we used the multiband VVV catalogs to determine 
extinction values by measuring the mean (/ - K s ) color of the 
RC stars and comparing it to that of a region of know extinction, 
named Baade's Window. A selection box in the CMD of each 
tile is chosen to include only Bulge red clump stars. The limits 
of this box vary for each line of sight depending on the amount 
of extinction and were therefore inspected by eye and modified 
accordingly. A Gaussian fit to (/ - K s ) color distribution is then 
used to determine the mean color of the red clump. For a 4' x 4' 
region in Baade's Window centered on coordinates 0:2000 = 18 : 
04 : 51.2 and 6 20 qo = -30 : 03 : 26.5, VVV photometry yields 
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Fig. 1. The mean value obtained from a Gaussian fit to a set of 
N data points randomly distributed following a Gaussian distri- 
bution centered on zero (solid horizontal line) and with a sigma 
of 0.10. The dashed vertical line shows the limit of 200 points 
adopted as the minimum number of stars in each subfield. 
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Fig. 2. VVV coverage of the Bulge in Galactic coordinates. Each 
small square shows the coverage by individual 1 ° x 1 .5° tile from 
the VVV. T ile numbering for the VVV survey is described in 
Saito et al.l (l2012h . The different color shading of the tiles de- 
pends on the spatial resolution (2', 4' and 6', respectively) used 
to derive reddening from the mean color of the red clump stars 
(see text for details). 



a (J - K s ) = 0.96 mag. Using the m easured extinction of this 
field of E(B -V) = 0.55 + 0.01 mag dSumill2004t Ikunder et alj 
2008; IZoccali etalJ [2008). we obtain a dereddened mean color 
of (J - K s )q = 0.68 for the RC stars in this region. We use this 
value as a reference for the intrinsic mean color of the RC, and 
its difference with respect to the measured RC color in any other 
region will be an indicator of extinction, meaning E( J - K s ). 



3. The complete high resolution extinction map 

We computed the mean RC (J - K s ) color in small subfields 
towards the complete Bulge region of VVV. We adapt the res- 



olution to provide a sufficient number of RC stars (> 200) to 
populate the color distribution with enough statistics, and at the 
same time to minimize the effect of differential extinction as dis- 
cussed in Paper I. The minimum number of 200 stars for each 
subfield was determined following the result shown in Fig. Q] 
The vertical axis of Fig. Q] shows the mean value obtained from 
a Gaussian fit to a set of N random values that follow a Gaussian 
distribution centered on zero and with a sigma of 0. 10. The ver- 
tical dashed line shows the limit of 200 stars. The dispersion on 
the mean value from the Gaussian fits increases significantly for 
a number of data points below this limit. The resulting resolution 
of the complete map depends on latitude, as shown in Fig. [2] For 
the latitude range -3.5° % b % +5° the density of Bulge stars 
allowed the highest resolution of 2' to be used. As the density of 
Bulge stars drops at higher distances from the plane, the resolu- 
tion had to be reduced to 4' for -7° % b % -3.5° and to 6' for 
-10° % b % -7°. The mean measured RC color is compared to 
the intrinsic value measured in Baade's Window to derive the 
reddening value E(J - K s ). From E(J - K s ), the Ak, extinc- 
tion values are obtained by applying a specific extinction law. 
Figure [3] shows the co mplete Ak, ex t inctio n map for the Bulge 
VVV region based on Carde Hiet al.1 ([1989) extinction law. The 
extinction variations in the inner 2° from the Galactic plane are 
better appreciated in the upper panel of Fig. |6]where the highest 
values of extinction reach Ak, =3.5. 

In a small number of regions the differential reddening is too 
strong even over 2' resulting in large errors for reddening mea- 
surements, or even in extreme cases preventing the identification 
of RC feature. Some of these regions are visible as small white 
dots/squares in the upper panel of Fig. [6] They are less obvious 
in Fig. [3] because of the resampling in the inner most regions in 
the presentation of the figure. All these regions are clearly iden- 
tified as having large errors in our extinction database. 

Our extinction maps show clearly the small scale variations, 
as produced by the strong dust features traced, in particular at 
low latitudes. These features are now seen in great detail thanks 
to the high resolution and large coverage of our study. 

3.1. A note about the extinction law 

At present, there is no real consensus on which is the correct ex- 
tinction law to be used for studies towards the inner Galactic 
bulge. Reddening maps are related to a given extinction law 
which can vary between different studies. Recent results suggest 
that extinction properties might vary in different locations and 
possibly depend on the amount of extinction and du st proper- 
ties. In particular. INishivama et "all (|2006[ I2008I I2009I) presented 
a detailed study of the interstellar extinction law towards the in- 
ner regions of the Bulge (|/| < 2°, \b\ < 1.0°) concluding that this 
varies significantly depending on the line of sight. 

INishivama et al.1 J2006) also provided ratios of total to selec- 
tive extinction which differ from a Ry ~ 3.1 standard extinction 
law as obtained towards other regions. 

The problem of the correct extinction law towards the 
Galactic bulge is not adressed here, but will be part of another 
article (Chen et al. 2012, in prep) based on GLIMPSE and VVV 
datasets. However, in order to obtain extinction values from 
the reddening map presented here, E(J - K s ) values need to 
be transformed into an absolute extinction through a given ex- 
tinction law. Here, we decided to adopt two different extinction 
laws. The firs t is the more commonly used extinction law from 
ICardelli et all (119891) . Adopting 1.240, 1.664, 2.164 micron as 
the effective wavelengths of a K2 giant (typical of red clump 
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Fig. 3. Extinction map of the Galactic Bulge for the comple te region covered by the VVV survey. The Ak s values are computed 
from E(J -K s ) measurements assuming Ca rdelli et al.l (1 19891) extinction law for all tiles. At Ak s values larger than 1 .5 mag the color 
scale saturates. The details of the extinction variation in the inner highly extinct regions, where Ak, reaches up to 3.5 mag, are better 
seen in the upper panel of Fig. [6] 



stars) for 2MASS fro mllndebetouw et alJ(l2005l) 
corresponding to the ICardelli et al.l (1 19891) extini 

A j = 1.692E(J-K S ) 
A H = L054E(J-K S ) 
A Ks = 0.689£(7 - K s ) 



, the coefficients 
extinction law are: 

(1) 

(2) 
(3) 



The second is the extinction law from Nishivam a et al.l (12009): 
Aj = l.526E(J-K s ) (4) 



A H = 0.855£(7 - K s ) 
A Ks = 0.528£(7 - K s ) 



(5) 
(6) 



which has been determined, for the 2MASS photometric system, 
in the inner regions of the Bulge closer to the Galactic plane. 

In order to test which extinction should be prefered, we per- 
formed the following exercise. We measured the colors (J - K s ) 
and (J - H) in our 4' control field in Baade's window, where 
we find mean values of (J - K s ) = 0.96 and (J - H) = 0.73. 
We then measure these colors for a set of regions with differ- 
ent latitudes and compare them to those in our control field. 
The E(J - K s ) and E(J - H) color difference for each line of 




3.0 



Fig. 4. (J - K s ) and (J - H) color difference between our con- 
trol field and those measured in the subfields of tiles b317, b303 
and b275. The size of the subfields corresponds to the same 
resolution described in Fig. 2. The blue dashed line shows the 
relation E(J - H) — 0.638E(J - K s ) c orresponding to the ex- 
tinction law from ICardelli et al.l (1 19891) and red solid l ine to 
E(J -H) = 0.67 \E(J - K s ) fromlNishiva ma et all (120091) 
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Fig. 5. Diffe rence between t he Ak , values obtained in this work 
and those of lSchlegel etall dl998), as a function of Galactic lat- 
itude, for 1000 randomly distributed 30' x 30' regions between 
(-8° < b < -1°). 



sight are shown in Fig. |4] The relation between these color dif- 
ferences can be obtained from the corresponding coefficients of 
each extinction l aw. Our measurements are in good agreement 
to those with the Nishiyama et al.l d2009l) extinction law, where 
E(J - H) — 0.671 E(J - K s ). However, since the implications 
discussed in Sect. 5 do not depend strongly on the adopted ex- 
tinction law, we proc eeded our analysis fol lowing the more com- 
monly used law from lCardelli et al" (Il989l) . 



3.2. BEAM calculator 

As described in the following, our high resolution extinction 
map of the Bulge provides a valuable tool for future studies. In 
order to make it available to the community, we developed a web 
based tool named BEAM (Bulge Extinction And Metallicity) 
calculatoiQ. This tool, provides the user with an easy access to 
our measurements, both for the extinction values presented in 
this article, as well as for the photometric metallicities obtained 
in Paper I. The complete metallicity map and its analysis will 
be subject of the forthcoming Paper III (Gonzalez et al. 2012, in 
preparation). 

The tool allows the user to input a set of coordinates and the 
corresponding size of the field of interest, retrieving the mean 
extinction {Ak,) towards these set of coordinates. The exti nc- 
tion values are comput ed using either ICardelli et al.l (1 19891) or 
Nishivam a et al.1 (|2009) extinction law. However, the E(J - K s ) 
values are always provided as well, so that the user can adopt 
any other extinction law of preference. 



4. Comparison with existing extinction maps of the 
Bulge 

Other studies have addressed the problem of extinction towards 
the Bulge. Reddening maps providing a similar Bulge coverage 
as the one pres ented here are t he 3D maps of Marshall et al.l 
d2006l) and the ISchlegel et al] d!998l) maps. As described in 
the introduction, it is well known that Schlegel maps are not 
trustable at regions closer to the galactic plane (\b\ < 5°) due 
to contaminating sources and uncertainties in the dust temper- 
atures. In order to quatify these effects, we compared our ex- 



tinction values to a set of lines of sight at different latitudes. 
We retrieved Ay values from the Schlegel maps for 1000 re- 
gions of 30' x 30' size, randomly distributed across latitudes of 
-8° < b < -1° for 0° < / < 2°. Although the relative dif- 
ference between Ag t values from Schlegel maps and ours will 
not depend on the adopte d extinction law, we c onverted their Ay 
values into Ak, following Cardell iet al.l (1 19891) law and then re- 
trieved our own Ak, values from the BEA M calculator for the 
same target fields, again selecting the ICardelli et all (1 1989b ex- 
tinction law. Figure [5] shows the comparison be tween our ex- 
tinction results and those of ISchleg el et al. d 19981) as a function 
of latitude. The difference its already visible at b < -6° and be- 
comes much larger for latitudes closer to the plane. Interestingly, 
we obs erve a negative difference between our Ak s values and 
those of Schlegeletal] dl998i) for b < -2° while it becomes pos- 
itive at b > -2°. The source of this sign variation is not obvious 
given the uncertainties in the Schlegel maps at those latitudes. 

Mars hall et al.l (120061) provide an extinction map for the com- 
plete Bulge, obtained by comparing the 2MASS photometry to 
the Besancon model, including the model dependent extinction 
as a function of distance in each line of sight. Unfortunately, 
this technique only allows to calculate the maps with a lower 
resolution (15'), which is insufficient to resolve the small scale 
extinction variations in the inner Bulge regions. 

Studies of lSchultheis et all (I1999L l2009h provided very high 
resolution maps (2'), obtained using DENIS and Spitzer data 
respectively for the inner \b\ < 2°. Since in these regions the 
small scale extinction variations become particularly strong, 
these datasets are ideal to compare with our extinction maps and 
study possible systematics be tween the methods. The Spitzer 
map of Schultheis et al. (2009) covers specifically the ~ 1° of the 
very central regions of the Galaxy. This region is also included 
in the DENIS map and the extinction derived from DENIS and 
Spitzer observations agree very well. 

In order t o pro duce the DENIS extinction maps, 
SchultheisetaD (Q999) used a method based on compari- 
son o f bulge RGB colors i n 2' subfields with model isochrones 
from iBertelli et al.l (Il994l) . They adopted a Solar metallicity 
population at a distance of 8 kpc with an age of 10 Gyr in 
the model isochrones. Same pro cess was adopte d to pr oduce 
the Spitzer extinction maps from ISchultheis et al.l (|2009). This 
method has the benefit of using the bright Bulge stars, thus 
minimizing problems arising from photometric errors and 
(in-)completeness. However, the colors of RGB stars are much 
more strongly affected by population effects, in particular 
metallicity and to less extent age, when compared to RC stars in 
our method. In fact, the color of RGB stars is frequently used to 
derive the metallicity distribution function of stellar populations 
when spectroscopy of individual stars is not available. 

The influence of the population effects on the RC 
magnit udes and color has be e n investigated both th eoret- 
ically dGirardi & Salarisl 120011: ISalaris & Girardil l2002|l_ and 



Pietrzyriski et al. 



observationally ( Alvesl I2000t iGro cholski & Saraiedini 



2002; 



20031: lLanev et al.1 l2012l) . In both cases. RC 



http://mill.astro.puc.cl/BEAM/calculator.php 



stars appear as reliable standard distance indicator, considering 
the possible variations in the general Bulge population proper- 
ties. If we consider a variation of the mean Bulge metallicity to 
go from Solar down to [Fe/H] = -0.4, as obtained f r om th e 
spectroscopic metallicity distributions of IZoccali et al.l J2008), 
we expect an error in the mean RC color of ~ 0.0 8 mag, equiva- 
lent to a variation of Ak, = 0.06 mag assuming the Cardell i et al 
(1 19891) extinction law and Ak, — 0.04 mag for Nishivama et al 
(2009) law. 
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Fig. 6. The upper pan el shows the inner ~ 4° region around the Galactic plane of our VVV extinction map. Ak, values are based on 
ICardelli et al.l d!989l) extinction law. Part of this region was also covered by the DENIS survey, used in Schultheis et al.l (1 19991) to 
build an extinction map shown in the lower panel. 



The extinction maps of lSchultheis et al.l dl999h . based on the 
magnitudes of RGB and AGB stars, were produced at latitudes 
b < 2° wher e no or only very weak metallicity gradient seems 
to be present (iRich et al.ll2007l) . Thus, extinctions obtained from 
isochrone comparis on to red giant stars sh ould be reliable for the 
regions studied by Schul theis et al.l d 1999b with errors dominated 
by the distance uncertainty of the RGB population towards each 
line of sight. However, this method would have a considerable 
problem when used to construct a global extinction map, such 
us the one presented here. 

Figure [6] s hows our extinction map and that of 
ISchultheis etail (fl 999) usin g DE NIS photometry. A v val- 
ues from ISch ultheis et all (11999) were converted to A% 
following the extinction law of ICardelli et al. (1989) for com- 
parison purposes. Their similarity is remarkable, with both 
maps tracing the same dust concentrations in the inner regions. 
We can also directly compare the ext inction values re t rieved 
fro m the BEAM calculato r to those of lSchultheis et al.l |l999) 
and lSchultheis et aTl d2009l) . This comparison is shown in Fig. [7] 
Ak s values of both methods are in good agreement with a scatter 
of the order of 0.2 mag, which is expected from the different 
methods and photometric errors. The scatter becomes larger, 
up to 0.6 mag, in the regions of higher extinction (Ak, > 2), 
where for a significant nu mber of fines of sigh t we find a 
lower extinction value than ISchultheis et ail (120091) . A possible 
explanation for this difference is that, in these highly reddened 
regions, the extinction is enough to bring RC stars to magnitudes 
fainter than K s ~ 17 where the completness of VVV is of the 
order of 60% (ISaito et al.ll2012l) and therefore our reddening 
values could be understimated. However, we emphasize that 



this difference is observed in the inner 1° x 1° of the Bulge 
where extinction determinations are certainly more complicated. 
Moreover, the overall agreement is remarkable considering the 
different techniques, data sets and stellar populations (RC vs 
RGB) of our studies. 

5. Implications for Bulge studies 

Extinction corrections are a fundamental parameter for essen- 
tially all studies of stellar populations, as well as for distances 
and the determination of structural parameters. In this section 
we highlight the advantages provided by our new high resolu- 
tion and large area extinction map for the studies of color mag- 
nitude diagrams, for determination of stellar parameters in spec- 
troscopic abundance studies and for distance and structure stud- 
ies of the Bulge. 

5.1. The Color-Magnitude diagram of the inner Bulge 

The VVV color magnitude diagrams are de-reddened using our 
new extinction map of the Bulge. In Figure [8] we compare 
the observed, extinction corrected, inner Bulge CMD with the 
predictions from the population synthesis model of the Milky 
Way from the Besangon group including the thin disk, thick 
disk and bulge pop ulations. T his model is described in detail in 
iRobin et al] d2012b . It uses the lMarshall et all d2006l) extinction 
map and when compared with the observations such as 2MASS 
star counts it can be used to study the structure of the Milky 
Way. As already mentioned above, in the innermost regions of 
the Bulge, 2MASS is severely incomplete and affected by strong 
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Fig. 7. Upper panel shows the c omparison of the Ay. values ob- 
tained in this work and those of Schu ltheis et al.l (I1999I) for the 
common inner Bulge region (-2° < b < +2°). Lower panel 
shows the comp arison with the SPITZER extinction map of 
ISchulfheis etaiT d2009l) for the central 1° of the Galaxy. Both 
plots are shown as density contours due to huge number of data 
points. 
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Fig. 8. Left panel shows HESS diagram of a 40' region in the 
inner Bulge (7, b) — (-1 .0, - 1 .0). The right panel shows the dia- 
gram for the same region as obtained from the Besancon model. 



crowding and blending, thus limiting the predicting power of the 
model in the inner regions. With the superior spatial resolution 
of VVV survey (median seeing is better than 0.9" in 7, H and K s 
images) and its deeper photometry, it is now possible to make 
direct comparisons with the simulations and thus validate some 
of the assumptions adopted in the model. 



Figure[8]shows a Hess diagram of a 40' x40' region centered 
at Galactic coordinates (l,b) = (1.0,-1.0) (VVV tile b320) in 
the left panel. In the right panel is the CMD obtained from the 
Besancon model for the same region. A remarkable similarity is 
observed between our dereddened data and the model. The ma- 
jor difference is for the blue plume of stars at (J - K s )o 0.2 

in the left diagram, which is much stronger, brighter, and bluer 
in the VVV CMD, than the blue sequence at (J - K S ) Q ~ 0.1 
in the model. This part of the CMD contains the disk stars, 
which are found all along the line of sight towards the Bulge, 
and therefore they have a range of reddening values, some hav- 
ing significantly lower reddening than derived for the Bulge red 
clump stars. Therefore after applying the extinction correction 
appropriate for the Bulge stars the colors of the disk stars can- 
not be correctly reproduced in our CMDs, but the colors and 
magnitudes of the major Bulge features, the RC and the RGB 
are in excellent agreement with the model. This provides evi- 
dence that VVV dereddened color-magnitude diagrams can be 
used to investigate the stellar populations properties, such as age 
and metallicity distributions of stars even in the inner Bulge re- 
gions. Additionally, the good agreement between the model and 
the observations indicates that we could use the Besancon model 
in comparison to observed VVV dat a, in order to obtain 3 D ex- 
tinction maps such as the one from Marsh all et al.l d2006). The 
high spatial resolution and photometric depth of VVV CMDs of- 
fer the possibility to derive the 3D extinction map of the Bulge 
with a resolution of up to 6' for the inner regions (Chen et al. 
2012, in preparation). Although this resolution, limited by the 
number of stars needed for the comparison, is lower than that of 
the map presented here, it will add the important 3D informa- 
tion, allowing more accurate characterization of the structure in 
the inner Bulge. 

5.2. Photometric effective temperatures 

Low and intermediate resolution spectroscopic studies, where 
determination of stellar parameters is based on photometric 
properties of the targets, depend on the adopted extinction cor- 
rections. In particular, the critical parameter is the effective tem- 
perature which is derived from the color of the stars. Current 
studies of the Bulge use multi-object fibre spectrographs with 
large field of view in order to observe simultaneously large num- 
ber of targets. The field of view can range from 25' radius for 
the FLAMES spectrograph on the VLT dPasquini et al.l 120031) 
to 2° field of AA Omega multifibre spectrograph at the AAT 
dSharp et al.ll2006l) . Assuming that a single value of extinction 
is adopted for a single pointing, using our 2' resolution extinc- 
tion map we can estimate the differential extinction in a typical 
region of say 30' and derive the effect on the corresponding ef- 
fective temperatures determined from colors of stars. 

The relative effect of differential reddening will not depend 
significantly on the adopted calibration, and in our experiment 
we use the photometric calibration from Ramirez & Melendez 
(2005) for giant stars. A typical Bulge K giant star will have 
Solar metallicity with an effective t emperature of around 4500 K . 
According to the calibration from lRamfrez &Melendezl(l2005h . 
this star will have an intrinsic color of (V - K s )q = 2.625. We 
now consider a 30' region located at a latitude of b = -4° along 
the minor axis (/ = 0°) and retrieve the extinction values using 
our BEAM calculator for a set of 200 randomly distributed po- 
sitions within this region. These individual extinctions will dif- 
fer from the mean extinction for the whole region depending on 
the amount of the differential extinction. If we add this differ- 
ence to the intrinsic color of the star and use the new color to 
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Fig. 9. The effect of differential reddening towards the Bulge for 
the determination of photometric effective temperatures. We plot 
the variations in the Ak s values for a set of 200 random positions 
in a 30' region towards b = -4° (black squares) and b = -2° (red 
circles), as obtained from the BEAM calculator, and compare 
them to variations in the photometric effective temperature. 



calculate the corresponding effective temperature from the pho- 
tometric calibration, we obtain an estimation of the effect of dif- 
ferential reddening when adopting a single reddening value for 
the whole region. The difference between the original effective 
temperature and the one obtained after including the effect of 
differential extinction is shown in Figure [9] as a function of the 
variation in Ak s . The amount of extinction Ak s varies by 0.10 
mag in the 30' region, which produces an error of up to 300 K in 
effective temperature. This effect becomes even more important 
in regions of higher extinction than Baade's Window as shown 
for the case at latitude b = -2° (red circles in Fig|9j». Here a 0.30 
mag extinction variation is observed within the 30' region which 
implies errors larger than 600 K in the derived photometric ef- 
fective temperature. 

Clearly all studies relying on photometric determination of 
the effective temperature need to take these extinction variations 
into account and correct for them. Otherwise the strong differen- 
tial reddening in the inner Galaxy may result in large errors for 
the abundance measurements. While the dependence of metal- 
licity error on temperature is not necessarily linear, and will de- 
pend on the spectral type of the target, for a reference we recall 
that already an error of 150 K in effective temperature led to a 
measu rement error in fFe/ Hl of 0. 1 dex for K giant stars in the 
Bulge dZoccali et alJl2008h . 

5.3. Distances and the Bulge structure 



As shown in iGonzalez et al.l d201 Id) dereddened magnitudes of 
the VVV survey can be used to build the Bulge luminosity func- 
tion towards different lines of sight and then measuring the mean 
magnitude of the RC to derive mean distance to different regions 
in the Bulge. This technique has been used in se veral studies 
to stu dy the orientation angle of the Galactic bar (IStanek et al.1 
119971: iRattenburv et all 120071: ICabrera-Lavers et al.l | 2007b and 



the X -shape morphology of the outer Bulge regions (Saito et al. 

12011b 

In IGonzalez et alJ (1201 lal) we studied the bar orientation 
in the inner regions (b = ±1°) by measuring the mean RC 
magnitudes as a function of longitude. Using the homoge- 
neous extinction map presented here, the study can now be ex- 




Fig. 10. Position of the Galactic bar with respect to the Sun as 
measured with the RC method from the the dereddened K s mag- 
nitudes using our extinction map (red filled circles) at a lati- 
tude of b = -5°. Also shown, are the measurements from the 
Galact ic model as presented in iGerhard & Martinez-Valpuestal 
(1201 2|) at the same latitude, corresponding to a bar with an an- 
gle of 25° with respect to the Sun-Galactic center line of sight 
(black empty squares). Dotted lines mark the lines of sight for 
longitudes 1= +5° and / = ±10° 



tended to other regions of the Bulge. Ideally, these measure- 
m ents should be compared to a ga lactic model such as the one 
of Martinez- Valpuesta et al. (2006). Interestingly, although this 
model has not been fin e-tuned in any way to reproduce the 
Milky way properties, in IGerhard & Martinez-Valpuestal (1201 2|) 
the model was shown to nicely reproduce the observed struc- 
tural properties of the Bulge, including the chan ge in the bar in- 
clinati on in the inner regions, as o btained by iNishivama et al.l 
(2005|) and IGonzalez et al.l d201 lal) . This model also provides 
predictions for higher latitudes which can now be validated 
using our extinction-free magnitudes. Figure [10] shows posi- 
tion of the bar with respect to the Sun as measured from the 
mean dereddened K s VVV ma g ntiude s and from the model of 
IGerhard & Martinez-Valpuestal (120121) at higher Galactic lati- 
tude b = -5°. The agreement between our observations and 
the model is remarkable, confirming that the Milky Way bar is 
inclined at 25 deg with respect to the Sun-Galactic center line 
of sight. The correction for differential extinction from our new 
map allows to extend this comparison of the model and observa- 
tions to the complete Bulge, including the characterization of the 
stellar density profile, the observed X-shape, and other structural 
indicators. 



6. Conclusions 

We have used the VVV data to obtain a high-resolution extinc- 
tion map of the Galactic bulge covering ~ 315 sq. deg area: 
-10.3° < b < +5.1° and -10.0° < I < +10.4°. We used the 
color of RC stars to trace the effect of extinction in sub fields of 
2' x 2' for -3.5° £ b < +5°, 4' x 4' for -7° < b < -3.5° and 
6' x 6' for -10° %b % -7°, providing a first reddening map sen- 
sitive to the small scale variations and covering homogeneously 
the whole Bulge. 

An excellent agreement is found when comparing the com- 
mon regions of our full VVV map with those of DENIS 
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(ISchultheis et al.1 1 19991) and Spitzer maps dSchultheis et al.1 
2009) which cover the very inner regions of the B ulge (\b\ < 2°). 
We al so have compared our maps with those of Schleg el et al.l 
(1998) at different latitudes and we confirm that Schlegel maps 
should not be used for latitudes \b\ < 5° towards the Milky Way 
bulge. 

The reddening map presented here becomes particularly im- 
portant when the global picture of the Bulge is investigated. We 
presented some of the direct consequences of using a higher res- 
olution map, when studying both the stellar populations and the 
structure of the Bulge. We show that a typical 30' region at 
b = -4°, for which normally a single extinction value is as- 
sumed, has a variation of up to 0.10 magnitudes only due to 
differential reddening. This effect becomes even more impor- 
tant in the inner regions where differential reddening produces 
variations of more than 0.35 magnitudes in a similar field size. 
This has a strong impact in stellar population properties obtained 
from the color of stars, such as the photometric effective temper- 
ature. We show that errors of more than 500 K can be caused by 
differential extinction in a 30' region if single extinction value 
is adopted. The current map allows to avoid such errors. Also, 
we show that the observed CMD of an inner bulge region at 
(l,b) = (1.0,-1.0) matches very well a CMD produced from 
the Besancon stellar population model where effects of extinc- 
tion are obviously not present. 

We have shown how a good treatment of extinction and the 
dataset such as the one presented here, provides the ideal tool 
for the inner Galactic structure study based on the luminos- 
ity function properties. We used the VVV data to measure the 
mean RC magnitude of the Bulge towards diff erent longitudes 
at b = -5°, following the method described in ([Gonzalez et al.l 
1201 lbl) . We find a remarkable agreem ent between our measure- 
ments and those from the model of Martinez- Valpuest a et al.l 
( 2006!) for a bar at 25° with respect to the Sun-Ga lactic center 
line of sight. iGerhard & Martinez-Valpuestal(l2012l) have shown 
how the model reproduces our observations in the inner Bulge 
(b = +1°) and we now show this is also true for larger latitudes. 
With our full extinction map it is possible to recover the global 
Bulge str uctural information in compariso n to models such as 
the one of Martinez-Valnuest a et al.l (120061) . 

Our extinction map is available to the community via a web 
based tool BEAM (Bulge Extinction And Metallicity) calcu- 
lator (http://mill.astro.puc.cl/BEAM/calculator.php). The com- 
plete photometric metallicity map (Gonzalez et al. 2012 in prep.) 
will also be available from the same site. These maps are an im- 
portant step towards the understanding of the general properties 
of the Galactic bulge and they are a particularly valuable tool for 
the forthcoming spectroscopic surveys of the inner Galaxy. 
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